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A review of some of the issues that have arisen over the years concerning the energy distribution 
among scales for magnetohydrodynamics (MHD) turbulence is given here. A variety of tools are 
employed to that effect, and a central role is played by taking into consideration the ideal (non- 
dissipative) invariants, namely the total energy, the magnetic helicity and the cross-correlations 
between the velocity and the magnetic field (concentrating on the three-dimensional case). These 
concepts, based mostly on theory, models and direct numerical simulations, are briefly put in the 
context of observations, in particular the solar wind, and some of the remaining open questions are 
delineated as well. New results on ideal MHD dynamics in three dimensions on equivalent grids of 
up to 6144 3 points using the Taylor-Green flow generalized to MHD are also mentioned. 

I. INTRODUCTION 
A. The context 

Astrophysical and geophysical flows are highly turbulent in general, be it only because they extend on a large ratio of 
interacting scales. Hence, one expects dissipative processes to play a secondary role as far as energy exchanges among 
scales are concerned, even though at the smallest scales where extreme events such as solar flares and coronal mass 
ejections can occur, heating prevails. Thus, it is thought that the overall dynamics results from nonlinear coupling 
between Fourier modes, waves, eddies and coherent structures in interaction. There have been many observations of 
turbulence in such flows. Our close environment such as the magnetosphere, the solar wind [221 11461 [157] and the Sun 
itself, as well as the interstellar medium (ISM) [3U] are particularly well documented for their turbulence properties; 
in the ISM, the high Mach number implies that interactions with the density and self-gravity play determining roles 
as well. Observations also indicate that magnetic fields are main agents of the acceleration of charged particles, 
contributing to shape auroral emissions and magnetospheric-ionospheric coupling, as in the case of Jupiter [128 . 

When such fields are present, they are often (but not always) in quasi-equipartition with either the gas pressure 
and/or the velocity. If plasma (kinetic) effects play a determining role at small scale, as observations show clearly [130J, 
the magnetohydrodynamics (MHD) regime - whereby the displacement current is neglected at low velocity compared 
to the speed of light - is important because it covers the dynamics of the large scales which hold the energy. Indeed, 
with an energy spectrum E{k) ~ fc _m , with 1 < m < 3 for energy and dissipation to be defined and interactions to 
be local, then most of the energy is contained at large scale described by the MHD limit. For that reason, however 
unrealistic this description becomes at small scales, the large-scale dynamics of complex magnetized flows can be 
described successfully by MHD. Beyond scaling laws, turbulent flows are known to develop characteristic structures: 
plasma sheet turbulence |153j and turbulence in a geomagnetic storm [129] were observed using remote sensing, such 
as with the CLUSTER configurations of four satellites, and such observations have yielded a wealth of detailed data 
on these complex flows, with for example coherent vortices |145j and vortex filaments in the magneto-sheath [3], Note 
also that the interplay between small-scale waves, such as Langmuir waves, and turbulent eddies contribute to forging 
the shapes of the wave-packets, as observed in the Solar Wind [62] . 

Even at high Reynolds number, the energy eventually gets dissipated if not in shocks in the supersonic case, then in 
intermittent and quasi-singular structures (vortex filaments, current sheets) that are as thin as the dynamics permits; 
roughly speaking, of the order of \/Rv, as for example for a one-dimensional Burgers shocks or a Harris current sheet 
described locally by a hyperbolic tangent profile, where Ry is the Reynolds number to be defined below. The precise 
shape of those dissipative structures, the fact that they roll-up or not (see [5TJ I108j for observations of such rolls in 
the Solar Wind), and hence presumably the rate at which energy is dissipated, in a finite time or not, does depend 
on interactions between large-scale and small-scale processes. There are special cases for which there is a proof for 
the lack of singularity in a finite time, such as neutral fluids in two space dimensions (2D), or fluids described by 
the Lagrangian-averaged approach in three space dimensions (3D) [251 (Ml 1MJ ; but in general, the available amount 
of energy in the system is transferred to the small scales where it is dissipated at the rate at which it is transferred 
(otherwise a run-off situation would occur). This is what has been observed rather unambiguously using direct 
numerical simulations, for neutral fluids [66l [68] , and in MHD both in 2D [Mj 1113) and in 3D [91] . 

Turbulent flows are known for their complex behavior, both in time (think of strange attractors) and in space. This 
is due to the nonlinearity of the advection term, and in magnetohydrodynamics (MHD), of the Lorentz force and 
Ohm's law. These quadratically nonlinear terms, under the assumption of incompressibility, lead, in Fourier space 
to a convolution which couples k, p and q modes together such that k = p + q is fulfilled; these are the so-called 
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triadic interactions which conserve energy and any other quadratic invariants Iq (see below) in each of their individual 
exchanges; hence one can truncate the system and still conserve the energy and the Iqs. 

In the absence of dissipation, nothing can arrest this coupling of modes, with the invariants being transferred both 
to large scales (as modeled by eddy noise or negative turbulent diffusivities for example, and as observed in inverse 
cascades in the forced dissipative case) and to small scales (as modeled by positive eddy viscosity and eddy resistivity, 
and leading to what is called a direct cascade). As far as jargon is concerned, let us reserve the vocable "cascade" for 
a transfer with a *constant* flux across scales in the so-called inertial range, as opposed to transfer which is simply 
referring to non-linear coupling with a possible scale-dependency of the flux. In this latter case, this may arise from 
Ekman friction for rotating fluids or from anisotropic Joule damping in the quasi-static limit of MHD at low magnetic 
Reynolds number [150, 151J. Also note that, in the ideal (non-dissipative) case, the cascades are flux-less. 

B. The Kolmogorov law 

The question we want to address now is: what is the distribution of energy among modes in a turbulent fluid? 
Anything? Or are there constraining rules? Phenomenology can come to the rescue: in the absence of any other 
time-scale but the advection time or non-linear time tjvl = i/u£ with ui the velocity at scale t, in the fluid case, 
there is a solution when one assumes, following Kolmogorov, that the energy spectrum depends only on the energy 
injection rate e = DEx/Dt and the wavenumber k, henceforth the so-called K41 law [70] : 

E(k) ~ e 2 / 3 fc~ 5 / 3 . 

This law is well verified in the laboratory and in observations of the atmosphere and the ocean, in certain parameter 
regimes, although there are corrections due to intermittency, i.e. the presence of strong small-scale localized structures 
in the form of vortex filaments. One can extend the K41 law to higher-order moments of the structure function of a 
given field u (say, the velocity field) computed on a distance r between two points, and assuming homogeneity: 

8u(r) = u(x + r) - u(x) , (1) 

with ujj = u • f the longitudinal structure function of u, the unit vector along distance r being defined as f = v/\r\. 
In the case of a self-similar field for classical (K41) turbulence, one has: 

(6u L (ry)~ r <> , ( p =p/3; 

this law is not verified by numerical data, nor by atmospheric data. Corrections to the self-similar scaling stem from 
intermittency, or the predominance at high-order of quasi-singular (dissipative) structures. 

One could also consider large-scale shear, in which case the same type of dimensional analysis gives E(k) ~ fc _1 
or, in MHD, one can think of the Alfven time associated with the propagation of waves along an imposed uniform 
magnetic field Bo, namely ta — £/Bq. The simplest solution could be to ignore these time scales and say that 
interactions will take place on the nonlinear time-scale as in the standard fluid case, i.e. the eddy turn-over time tjvl- 

Is this the case? 

This paper is devoted to a rapid review of where we stand now on this thorny problem. Thorny because it is difficult 
to measure spectra in astrophysics whereas it is easy to do so, at least in principle, with direct numerical simulations 
(DNS) under the assumption of isotropy. However, the Reynolds numbers achieved in DNS in three space dimensions 
are moderate, hence the extent of the range of scales where such laws can apply omitting the influence of the forcing 
and of the dissipation, is not very large. And this issue is thorny because, when there is more than one time-scale 
involved, dimensional analysis alone cannot give us the answer. 

Moreover, in the case of forced flows achieving a statistically steady state, the time-scale associated with the forcing 
itself (white noise, red-noise with some coherence time, or constant) may play a role as well. So could the existence 
or not of a uniform imposed magnetic field, as well as the aspect ratio of the computational box and the presence of 
boundary layers. Also, it is found in numerous studies of two-dimensional fluid turbulence [MJ 11061 1149] that when 
adding friction at large scale, this can modify both the aspect of the coherent structures that form in the inverse 
cascade and the scaling law for the distribution of energy (such structures can be detected using wavelet algorithms, 
see [78] and references therein). Similarly, there are observations, both in the solar wind |112) and for solar flares [T], 
that indicate variations with time/space of spectral exponents for MHD turbulence. 

Finally, one needs to stress that the value achieved by such a spectral index in MHD turbulence is not simply a 
matter of a theoretical nature or numerology. It can also have consequences, such as on the resulting heating rate of 
coronal events as a function of the axial field |127j . 
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C. The equations 

We give for reference the dynamical equations for MHD. They read, for an incompressible fluid with v and b 
respectively the velocity and magnetic fields in Alfvenic units, assuming a uniform density equal to unity: 

— +v Vv = -VP+j x b + i/V 2 v, (2) 
ot 

— = V x (v x b) +y?V 2 b . (3) 
at 

Note that b is, dimensionally, a velocity as well, the Alfven velocity. V is the total pressure, V ■ v = V • b = 0, and 
v and 77 are respectively the kinematic viscosity and magnetic diffusivity as stated before. With v = 0, i] — 0, the 
energy Et, the cross helicity He and the magnetic helicity Hm, defined as 

E T = E V + E m = (v 2 + b 2 ) /2 , H c = (v b) /2 , H M = (A • b) /2 , 

where b = V x A with A the magnetic potential, are conserved (factors of 1/2 in the helicities are here for convenience, 
to agree with definitions given in |42|1. Note that in the strict two dimensional case (2D, v z = 0, b z = 0), magnetic 
helicity is identically zero and this invariant is replaced by (A 2 } (and higher order moments associated with so-called 
Casimirs, which are not preserved by the Fourier truncation, though, beyond second order). 

The kinetic energy spectrum is defined as the Fourier transform of the velocity two-point correlation function. Once 
homogeneity, isotropy and incompressibility have been taken into account, only two defining functions remain: Ey(k) 
is proportional to the kinetic energy, with J Ev(k)dk = Ey = |(t' 2 ), and the kinetic helicity, Hy(k), stems from the 
anti-symmetric part of the velocity gradient tensor; helicity is a pseudo-scalar, with J Hv(k)dk = |(v • w); Hy is an 
invariant in the absence of viscosity in the pure fluid case but not in MHD. 

One can also introduce the Elsasser variables = v ± b; defining D±/Dt = dt + 1 ■ V, one can then obtain a 
more symmetric form of Eqs. ^ and ^ and write: 

D z rZ ± V + 71— 9 _(_ V~ 77 9 _. ... 

-gp =-^+V vV +V vV - (4) 

There are no self-advected non-linearities z + z + or z~z~, hence the common conjecture that in some sense MHD 
turbulence is weaker than its fluid counterpart, the effect being associated with Alfven waves (corresponding to 
z ± =0). Finally, the kinetic and magnetic Reynolds numbers are defined as 

R v = U Lq/v , R M = U Q L /r), 

where Uq, Lq are the characteristic velocity and length scale; one could define similarly 

R± = Z^L Q /v ± , 

with = Uq ± Bq and 2v± = v ±r), with Bq being the magnitude of b, a characteristic large-scale magnetic field. 
We shall assume v — r\ in the following, hence Ry = Rm and the magnetic Prandtl number Pm — v/rj — 1 unless 
otherwise stated; however Zq can be quite different, depending on the amount of cross-correlation in the flow. Of 
course, when v = r], hence v_ = 0, the ± Reynolds number defined above reduce to one, R+. 

In the next section, we analyze the case of ideal (non-dissipative) dynamics in three space dimensions, and examine 
energy spectra in §III| The role of the exact flux laws that can be written in MHD is reviewed in { IV and an example 



of lack of universality in MHD is dealt with in §Vl The next two sections examine the role played by the degree of 



alignment between the velocity and the magnetic field (j |VI[ ), and by magnetic helicity (SVII). In the penultimate 
section, one deals briefly with how models of turbulence can help, and finally §IX| is the conclusion. Note that some 
of the concepts to be delineated below are well known, see e.g. [13] (see also |28 [ I123ffl26j . and references therein), 
and the implications for astrophysical turbulence have been reviewed for example in [351 11461 1157j . 



II. WHAT CAN WE LEARN FROM STATISTICAL MECHANICS 



A long-standing concept is that the dynamics of turbulent flows is fully encompassed in the nonlinearities of the 
primitive equations, provided we are not in a weak turbulence regime where waves may be the dominant effect, at least 
at early times, and provided as well that the Reynolds number is very high, i.e. that, for the most part, dissipative 
effects are negligible. Note however that, through mode coupling, there is no way to arrest the cascade to small scales 
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and therefore there exists a scale, called the dissipative scale in, at which dissipation sets-in, no matter how small the 
viscosity is; it is found simply by equating the time associated with dissipation, Tdi sa = ^dI v j an d the relevant (say, 
nonlinear) time at that scale. This leads to what is called reconnection of field lines that are broken by dissipation (see 
e -g-> [7^1851 1581 II 131 ITT51 11341 H4"8] ). It gives rise, in astrophysics and plasma physics, to energetic phenomena such as 
coronal mass ejections and solar flares (see |104j for a large-scale reconnection event observed in the magnetosphere) . 

In such a framework, one can consider the dynamics of a truncated ensemble of modes under their nonlinear 
dynamics. It was shown by T.D. Lee |75j . for three-dimensional fluids and MHD, that one can write in fact statistical 
equilibria for such systems, taking into consideration the only constraint, in the absence of dissipation, of the energy 
conservation since it is preserved by the truncation (for the one-dimensional problem, see [331 E7J and references 
therein). In the simplest case, this leads to equipartition among Fourier modes, and thus to an energy spectrum in 
dimension d, E(k) ~ fc^" 1 , i.e. proportional to the number of modes in each (isotropic) Fourier shell. 

It was soon realized that, when there is more than one invariant (in fact, the generic case for fluids and in MHD as 
well), strange things can happen: namely, that when the relative importance of these constraints, due to invariants, 
differ, the resulting spectra can vary profoundly. On this basis, Kraichnan postulated an inverse cascade of energy 
toward large scales in 2D Navier-Stokes, assuming that adding viscosity to the system would allow for an energetic 
balance but would not alter the non-linear dynamics (see [72 ). It was shown recently for the Kolmogorov spectrum 
in 3D ideal fluids obtained at intermediate times and intermediate scales [33], that the small-scale k 2 equipartition 
spectrum acts as an effective eddy viscosity to these intermediates scales, even though the equipartition spectrum is 
progressing, slowly, toward the large scales. Even in the presence of waves, one can observe these equilibria although 
it takes a substantially longer time as the linear term grows in strength, for example in rotating turbulence as the 
Rossby number decreases [93] (note that the so-called "1/f" noise can also develop in turbulent flows for long times, 
showing the importance of memory effects, see e.g. for a recent study |38j and references therein). A similar conclusion 
of intermediate time-intermediate scale turbulent-like behavior in ideal fluids is reached in MHD in two dimensions 



[4"31 [731 [152] . although the power law index is not necessarily a K41 law (see ^ III ) . Presumably, the same happens in 
three dimensions as well (see also [58]). 

So, what are these statistical equilibria in 3D MHD? They were derived in [32]; with a, (3 and 7 the Lagrange 
multipliers associated with the Et, Hm and He invariants, namely 

aE T + pH M + 1 H C ; 

note that j3 does not have the same physical dimension as a and 7. These equilibria read: 

H J (k) = k 2 H M (k) = -^ i ^ ) ; H G (k) = ^Hj(k) , H v (k) = ^ Hj(k) , (5) 



E M (k) = -^H j{ k) = -g^ , M *)=(i-^i)r w .A.. 



where Hj = J k 2 HM{k)dk is the current helicity, and 



a>0, r 2 = l-^>0; V{k)= U-^I^J >0Vfe€ [k min ,k max ] ■ (7) 

Conditions on these multipliers and relations with the minimum wavenumber in the system are such that real- 
izability is ensured (see |42] for details). By realizability is meant positivity of the energy spectrum, or a Schwarz 
inequality relating energy and helicity spectra, namely H (k) < kE{k) with maximal helicity in the case of an equality, 
corresponding to alignment (parallel or anti-parallel) of say the velocity and the vorticity. 

All Fourier spectra are proportional, i.e. equivalent, to within constants based on the Lagrange multipliers i.e. on 
the relative importance of the three invariants imposed to the truncated system of modes at t = 0. This is true except 
for the kinetic energy, which has to recover its non-MHD formulation when b = 0. More importantly perhaps, Vfc, one 
has Ey(k) < Ej^ik), a point already noted in [143] , In fact, defining the relative modal energy and helicity Eu(k) 
and Hn(k), we have En(k) < 0, and flfl(fc) * HM{k) < 0, or more specifically: 

E R (k) = E v (k) - E M (k) = -E M (k)^^ <0 , (8) 



H R (k) = H v (k) - k 2 H M {k) = -T 2 k 2 H M {k) 



(9) 
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Examining a bit further these statistical equilibria (see also |141j ). one observes that there can be an inverse cascade 
of magnetic helicity for proper values of j3 with respect to the other two multipliers a, 7 (note that Hq^m are pseudo- 
scalar and that, unlike energy, they can have either ± sign). Moreover, equipartition (E^ = 0, Hr = 0) does not 
occur except for zero magnetic helicity (/3 = 0), or for maximal correlation in the initial state, as well as for very large 
systems (k — > 00). It should be noted here that a lack of equipartition is viewed consistently in the solar wind 83 , 
with most often a slight excess of magnetic energy. It can be interpreted as due to strong v-b correlations, arising from 
the source. However, since a similar trend is also seen in numerical simulations and in models of MHD turbulence, 
in each case the slight excess of Em can likely be attributed to ideal dynamics prevailing at high Reynolds number. 
Other remarks, perhaps a bit unexpected, are that: 

• In the absence of correlations (7 = 0, T 2 = 1), of course Hc{k) = 0, but so is the kinetic helicity, Hv(k) = 0, 
which can thus be viewed in ideal MHD as being due totally to the amount of alignment between the velocity and 
the magnetic field (or between the magnetic field and the magnetic potential). In that case, the kinetic energy 
spectrum has its fluid value, Ey{k) — Airk 2 /a (the magnetic helicity playing no role), whereas the magnetic 
energy is EM(k) — Ev(k)/T>(k); so that, for wavenumbers approaching |/3|/a, magnetic energy becomes very 
large, and so does magnetic helicity. This is at the origin of the inverse cascade excitation of large scales in 
MHD turbulence, as postulated in [32] ■ For 7 7^ 0, such condensation of magnetic excitation can happen as well. 

• In the opposite case of maximal correlation, equipartition is recovered (Eji(k) = 0, Hn(k) = 0), as expected. 

• Furthermore, note that, rewriting in part equations (5-7), we have: 

H c (k) = -J-E M {k) Vfc . (10) 



In other words, when 7 7^ 0, the correlation spectrum in ideal MHD is proportional to the magnetic energy 
spectrum, the latter being constrained by a Schwarz inequality to be larger than the magnetic helicity, namely 
Ei\i(k) > kH]\i(k). Thus, when magnetic helicity undergoes an inverse cascade to large scales, magnetic energy 
has to grow to large scales as well, and in the force-free (fully helical) case, it will have a spectrum determined 
by that of magnetic helicity; and thus the correlations will grow as well to the large scales. A negative transfer 
at large scale was observed in some early low-resolution numerical simulations of the MHD equations |122j ; this 
remark could also explain recent observations in the Solar Wind [140] . 



It should also be noted that the correlation spectrum docs not change sign at any scale, as evident from eq. ( 10 1. 
Such a change of sign can thus only occur through forcing mechanisms or else dynamically at the dissipative 
scale, as it has been observed and modeled using phenomenological arguments in two-point statistical closures 
of turbulence [57] (see [50] for 2D DNS). Similarly, kinetic and magnetic helicity are of the same sign at all 
scales, as Eq. ^ indicates; this sign is determined by the initial conditions and the helicity spectra can only 
change sign beyond the dissipative scale (leaving aside here a discussion involving the case when the magnetic 
Prandtl number Pm — v/rj 7^ 1 for which more than one dissipative scale can occur). It is not clear whether 
this latter point concerning the behavior of 7?Aj-(fc) at the dissipative scale has yet been verified by data. 

Other cases of ideal MHD flows have been studied, and the invariants change as one either changes space 
dimension [82 , 98 , or adds either a uniform magnetic field B , or a uniform rotation Q ; in the case when 
Bo 7^ and £!o 7^ 0, a new invariant appears |136l 1137] , called the parallel helicity 

Hp = H c — a H M , 



when these two externally imposed agents are parallel, namely with f2o = ffBo, even though He and Hm are 
not invariant anymore in that case (see [138] for the corresponding two-dimensional geometry). 

Finally, ideal dynamics is of course well-suited for studying the possible development of singularities in a finite time 
in MHD, a problem that is open, as it is as well in the three-dimensional fluid case. Although this question is of 
a mathematical nature, one can help with direct numerical simulations that are accurate so that the invariants are 
well preserved. One such example of ideal 3D MHD dynamics is given in the following figures. In the computation 
described succinctly below, there is no dissipation, no forcing and no imposed uniform magnetic field; the magnetic 
helicity is identically zero, and so is the cross-correlation. Initial conditions for both the velocity and the magnetic 
field have the same four-fold symmetries as that of the Taylor-Green vortex [17] but extended to MHD [76l [77] . 
Such a vortex corresponds to the flow between two counter-rotating cylinders as encountered in many experimental 
configurations. One can compute the evolution of such flows until the dynamics reaches the mesh size Ax — 2tt/N 
where N is the number of points per dimension. After that time, the partial differential equations are no longer 
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FIG. 1: Temporal evolution of the logarithmic decre- 
ment S(t) (see text, eq. (Ill) for the ideal MHD 



case with, as initial conditions, the so-called I-flow 
(see [76 ). The horizontal dashed line indicates the 
mesh size Aa;, on a run using the equivalent of 3072 3 
points, implementing the symmetries of the Taylor- 
Green flow generalized to MHD. 



FIG. 2: Temporal evolution of the energy spectral 
index n{t) (see text, eq. ( 11 1) for the same run as in 
Fig. [l] The computation of the MHD equations is 
reliable until 8 ~ Ax, indicated by the dash line in 
Fig. [I] or t w 2.5, time after which one enters the 
regime of ideal dynamics eventually described by the 
statistical spectra given in Eqs. (5-7). 



resolved and one deals with the dynamics of a truncated number of Fourier modes which will evolve progressively 
toward the ideal spectra given in equations (5-7). 

When fitting the temporal evolution of the energy spectra as: 

E(k, t) ~ C{t)k~ n ^ eyip [ ~ 2Sit)k] , (11) 

where S(t) is called the logarithmic decrement and n(t) is the inertial index, one is assured of the regularity of the 
MHD equations as long as 8 ^ 0. In Figs. [T] and [2] are displayed the temporal evolution of 8 and of n for an ideal 
MHD run on an equivalent grid of 3072 3 points; the horizontal dash line represents the mesh size which is reached for 
t « 2.5; just before that time, a sharp acceleration of the evolution of 8 is observed, as in the computation reported in 
|76j where it was interpreted as the formation of a rotational discontinuity between two approaching current sheets. 
After that time, the computation of the PDEs is not reliable any longer. A higher-resolution run on an equivalent grid 
of 6144 3 points has allowed for the further study of this evolution [18] . The dominant structure that emerges consists 
in the collision and subsequent coupled evolution of two neighbouring current sheets. At the last reliable time of the 
computation (i.e. when the logarithmic decrement is comparable to the grid size), magnetic field lines embedded in 
this dual system of current sheets are all seen converging to the point of maximum current [18] . 

Note that the spectral index of the total energy spectrum seems to reach values below n — 3 (see [7511152) for recent 
studies of the corresponding ideal problem in two dimensions). Finally, in Fig. [3| are given the Fourier spectra of 
the total energy, with different successive interval of times displayed in different colors, and with a change in color 
every AT ~ 0.4. At the latest time displayed here, the spectrum begins to show an accumulation at small scale, 
the premises to the equipartition spectrum at small scales following a k 2 law in the simplest (non-helical) case. The 
jagged aspect of the spectra at large scale is likely due to resonances between modes because of the symmetries of the 
initial conditions. 

However, the analysis of ideal dynamics in MHD docs not tell us what actually happens in the inertial range at 
intermediate scales for MHD turbulence in the presence of dissipation, nor does it put any constraint that we can think 
of in the inertial ranges themselves. So we now turn to analysis and phenomenology, using dimensional arguments 
based on characteristic time scales, in order to advance in our understanding of MHD dynamics. 



III. ENERGY SPECTRA 



A. Weak MHD turbulence 



Let us assume that Sb/B << 1 where 8v and 8b are velocity and magnetic field fluctuations in the presence of a 
strong uniform magnetic field Bp. Under the effect of such a strong imposed field, Alfven waves will propagate, with 
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FIG. 3: Temporal evolution of the total energy spectrum Er(k, i) for the same flow as in 
Fig. [I] and Fig. [2j in logarithmic coordinates. Colors shift from blue to green, red, yellow, 
black, purple and cyan with increments of roughly AT = 0.4. The first spectrum is at 
t = 0, and the last one is at t = 2.815, time at which pile-up of energy in the smallest scale 
is clearly visible; it would lead at later times to an energy spectrum ~ k 2 , corresponding 
to the non-helical ideal MHD. 



5v ~ 5b. Hence, the eddy turnover time is slow compared to the Alfven time based on Bq » Sb, and we have a small 
parameter in the problem, namely: 

/ Sb h 

ewr = ta/tnl = -5-7- j 

assuming the turbulence becomes quasi bi-dimensional, with £y and £± referring to directions parallel and perpen- 
dicular to Bo- This is a situation familiar in weak turbulence theory: the nonlinear interactions are weak except 
when the two counter-propagating waves are at resonance, and thus the statistical problem in terms of moments and 
cumulants of the fluctuating fields can be closed (see e.g. for review [103] and references therein). This was exploited 
in the case of MHD in pHl |4"5] . The procedure leads to integro-differential equations for the diverse spectra (kinetic 
and magnetic energy and helicity) and such equations can be analyzed for their flux-less solutions (corresponding to 
statistical mechanics of ideal truncated systems), and for their constant flux solutions, corresponding to the weakly 
turbulent regime. In the latter case, one finds 

£(fc ± ,fcn)~fc: 2 /(fcii) (12) 

(henceforth, the WT spectrum), where the perpendicular and parallel directions refer to the imposed magnetic field. 
Note that, at this order, weak MHD turbulence can be viewed as being degenerate, since there is no coupling between 
different planes parallel to the magnetic field axis. One advantage of the weak turbulence approach is that not only 
can one determine the power-law scaling but also the (generalized) Kolmogorov constant appearing in front of the 
scaling which is known as well; it is found to depend strongly on the degree of alignment between the velocity and 
the magnetic field [44] . 



Phenomenology can help us in recovering this spectrum in eq. (12 1, with "proper" dimensional analysis, by using 
ewr in a hand-waving argument. In the presence of a strong Bq, the turbulence is weak, hence the nonlinear transfer 
is weak; or perhaps, better said, it is delayed. By how much? Let us use the only parameter we have at our 
disposal, namely ewT- Hence let us conjecture, following Iroshnikov (1963) and Kraichnan (1965) [5SJ [7T] that the 



characteristic time of transfer to small scales ttr is longer than in the fluid case by 1/ewT, hence ttr = t^ l /ta- 
Dimensional analysis then gives, with t^l = £±/ue ± using the fact that most of the nonlinear transfer occurs in 
planes perpendicular to B , and ta = £\\/B Q following the anisotropic dispersion of Alfven waves, 

E(kx,k\\) ~ k- 2 k~ 1/2 . 

In the isotropic case (k± ~ hi), and after integration over one direction, one has E(k) ~ [e T 5o] 1/2 fc _3/2 , henceforth 
the Iroshnikov-Kraichnan (IK) spectrum [55] [7T], with er = DEt/DT the rate of dissipation of the total energy. 
Hence the IK spectrum is compatible with the theory of weak MHD turbulence. 



B. When the weak regime breaks down 



Are such spectra observed? This is where trouble begins. Of course, in the weak turbulence integro-differential 
system of equations, they are obtained analytically, and they are consequently observed when performing a numerical 
integration of the complex equations emanating from the weak turbulence development [JH US] . But the weak MHD 
theory is non-uniform in scale, because the two characteristic time scales vary in different ways with the size of eddies; 
hence, there exists a scale at which cwt ~ 1 and some sort of strong turbulence takes over, to be determined. This 
is why it is rather difficult to observe such a spectrum although there are indications that it has been observed in 
the magnetosphere of Jupiter [131] for which 8b/ Bq ~ 0.008 and one estimates cwt ~ 0.06. Saturn could also be a 
candidate for such a regime to be observed. Note that this situation of break-down of weak turbulence is well known 
in geophysical fluid dynamics, where the Ozmidov scale for stratified flows and the so-called Zeman scale for rotating 
flows have been defined in similar ways (as the scale for which the appropriate ewT = !)• It was shown recently 
using a computation on a grid of 3072 3 points of rotating turbulence that a Kolmogorov spectrum recovers beyond 
the Zeman scale; such a high resolution was needed in order to be able to resolve the large scales, the first (weak 
turbulence) inertial range, the K41 range and the dissipation range as well |94j . 

In the weak turbulence regime, Alfven waves are fast and equipartition obtains at all scales, except possibly at the 
largest scales where magnetic helicity would dominate. However, it is well-known that a defect of equipartition has 
been observed in Solar Wind data, with what is called an (inverse) Alfven ratio EM{k)/Ey{h) often found to be close 
to 2, similar to models and DNS of MHD turbulence. We show in Fig. [4] such a ratio at the peak of dissipation in 
several runs done on moderate grid resolution |141j : the spectra are seen to differ at large scale but are close to unity 
in the inertial range, as already found in 77J. This numerical result is reminiscent of Solar Wind observations of the 
Alfven wave ratio that also appears to be constant in the inertial range [35] . 

There are impressive observations in astrophysics which exhibit, combining different techniques, a Kolmogorov 
spectrum for the density (and hence, probably for the energy as well), on many orders of magnitude, from the 
astronomical unit to the size of the inter-galactic medium [6] . Of course, errors are large, and other spectra could 
be accommodated as well by the data. Note that one reason why such spectra might be observed is because even 
though this medium is in general supersonic, there are strong pockets of subsonic flows where K4I spectra can develop 
as indicated by numerical simulations, for example using modeling of the small scales. In the 
Deen known for a long time that the energy spectrum is often quite close to K41 [83] (see also the 



Solar Wind, it has 

reviews [H [US QUI), but not always by far [TT2] , 

One could argue that an anisotropic Kolmogorov law is the solution to MHD turbulence, once one takes into 
account that the main effect of a strong magnetic field is to render the turbulence quasi bi-dimensional. This is 
basically the solution advocated in [53]. The argument used is one of "critical balance" at all scales between Alfven 
wave propagation on a time ta(£) and non- linear transfer on a time tnl(£) at scale £, namely that at all scales, one 
has eGs(£) — ta/tnl — 1, leading to hi ~ fc 2 / 3 / B (see [15] for a recent review). This scaling is corroborated by some 
numerical simulations [HI] and not by others [771ll01lll41j . Moreover, note that one can extend this phenomenological 
argument to a constant ratio with scale, but not equal to unity, i.e. ecs2 — c W, I, in which case other solutions 
(such as the IK spectrum and the WT spectrum) can emerge as well [3S] . The analysis of data is further complicated 
by the fact that, for different angles between the imposed field Bo and the wavevector k, different spectra may occur, 
as documented in rotating flows [M] and in MHD as well [TT], with the slow modes (corresponding to zero-frequency 
standing waves) and three-dimensional modes evolving differently, the latter having a steeper spectrum. Henceforth, 
averaging over angles as is done to obtain isotropic spectra may not be justified (see also [T2])- 

Recent numerical simulations taking care in obtaining convergence as the resolution is increased conclude that the 
GS model is what MHD will tend to, as the Reynolds number is increased [§]. However, it should be noted that the 
runs with hyper-diffusivities (dissipation operator with orders higher than a Laplacian) create their strong bottleneck 
effect, and that from an examination of the data in [5] one may not yet be able to discard IK spectra in these instances. 



[S] (see also SVIII 
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FIG. 4: Ratio of spectra of magnetic to kinetic energy at peak of dissipation for several 
runs in three dimensions (see also [141] ). Note the differences at large scale, presumably 
responsible for different scaling laws, and the constancy in the inertial range, with a slight 
excess of modal magnetic energy as also observed in the Solar Wind 35 . 



C. Inter mittency 

Intermittency deals with rare but strong events in turbulent flows, as well as in other critical phenomena [20] , 
Because of this phenomenon, the data for the energy spectrum in MHD turbulence is perhaps still a bit ambiguous: 
after all, the K41 spectral index 5/3 is not very far from the IK solution of 3/2 in the isotropic case, and besides 
we expect these values to be different because of intermittency corrections due to the presence of strong localized 
structures, in the form of vortex filaments for three-dimensional fluids, and vortex and current sheets |113[ 1115] 1148] . 
that eventually roll-up in 3D [55] , as observed in the magnetosphere [BT1 1108] (although in that case the roll- up is at 
times interpreted as being due to Hall currents). 

When examining higher-order moments of structure functions, the ambiguity is completely lifted: MHD turbulence 
does differ rather strongly from fluid turbulence in its scaling laws, independently of isotropy. In two dimensions 
[118] , as well as in three dimensions [SHHHH], intermittency exponents are very different from the fluid case, and MHD 
is more intermittent (in the sense that there is a stronger departure from self-similar linear scaling, although in that 
latter 3D case, to a lesser extent than for 2D). Such intermittency can be described by somewhat ad-hoc models, 
following the lead of She and Leveque [135] , as done for example in [501 1114) , and a link to self-organized criticality 
has been made recently for 3D numerical simulations of MHD turbulence |147] . 

Similar results hold, moreover, in observations of the Solar Wind [23] and more recently in auroral absorption |142j 
as well as in X-ray emissions from the solar photosphere in active regions 11], with quite dramatic, order unity, changes 
in scaling exponents when going from weak flares (for which the intermittency is quite close to that of neutral fluids) 
to large flares in which case the exponents are closer to the MHD case in 2D, presumably indicative of the presence of a 
strong quasi-uniform field. Note that one also finds different intermittency exponents for supersonic MHD turbulence 
when varying for example the Mach number [105] . 

Is it thus conceivable to claim that only one solution will emerge for energy spectra at high Reynolds number, 
considering that the Reynolds number of astrophysical flows in which variable intermittency is observed, such as in 
the Sun, is quite high already? Can these differences only be ascribed to measurements errors (granted, they are 
large)? In particular in light of the fact that the absolute difference in anomalous scaling indices between strong and 
weaker solar flares is of order unity? Or is it the difference in physical parameters (compressibility, radiative transfer, 
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stratification, boundaries, geometry) that is the culprit? Perhaps not. So, what can we do? 



IV. EXACT LAWS 



There is one thing we know about MHD turbulence: there are exact laws that can be written, under hve simplifying 
hypotheses, and such laws can be of some help. Under the assumptions of incompressibility, homogeneity, stationarity, 
isotropy and high Reynolds numbers, one can write for the scaling of third-order structure functions defined before 
in eq. ([l]), for the Elsasser variables, z ± = v ± b |116[ 1117] : 

(SzZ(r) £«[^(r)] a ) = - J e + r ; (<5z+(r) £^r(r)] 2 ) = -*e~ r , (13) 
or, in terms of the velocity and magnetic field: 

( Sv L ^(Svi) 2 ) + {5v L Z t (5b t ) 2 ) - 2( 5b L ^S Vl Sbi ) - -~ e T r , (14) 



( 5b L Xi(8bi) 2 )-(5b L S,;(^) 2 ) + 2( 5v L EiSvi Sb t ) = 



4 c 
— : e° r 



(15) 



Such laws stem from the conservation of energy and cross- helicity. A detailed study of the relative importance of these 
terms has been performed in |156| for a set of numerical simulations. Similar laws exist for Navier-Stokes turbulence, 
and for the dynamics of passive tracers in fluids. Note that these exact laws have been extended to several cases, for 
example with anisotropy [48], or including the Hall term [47] , 

One essential remark is that equations (14) and (151 are coupled: the invariance of total energy (governed by 
the imposed flux of energy ex) and the invariance of cross-correlation, governed by its rate of injection ec chosen 
independently of ex, play dual symmetric and coupled roles; it is particularly visible when considering the Elsasser 
variables, and their rates e± = [ex ± e c]/2- Each rate can impose its own time-scale, say Tc for the energy and 
cross- helicity (or T± ) , and thus there is a priori the possibility of breaking the universality that is often assumed in 
the neutral fluid case. Note that the possible role played by the invariance of kinetic helicity for neutral fluids is 
discussed in [SS], and the case of magnetic helicity in MHD is analyzed separately in { VII (see also [119] ). 

Let us examine a bit further equations (14) and (15). There are three possibilities a priori: either Sv >> 6b (the 
kinematic case of the dynamo starting from a weak magnetic field, in which case we certainly expect a Kolmogorov 
law for the total enegy spectrum dominated by its kinetic part), or Sv << 5b (corresponding to a strongly magnetized 
force-free plasma), or of course Sv ~ Sb, which can be called the Alfvenic case. Such partition of phase space has been 
documented, both in 2D and 3D numerically |141[ 1 1431 1144) and is plausible. But what do the exact laws above tell 
us in each case? 

In the kinematic regime, one recovers approximately the scaling for standard fluid turbulence, with 
(SvL(r)'Ei5vi(r) 2 } = — |er (in dimension three), compatible with the perhaps more familiar expression in terms 
only of the longitudinal structure function (5vL(r)) 3 = — |er, with a Kolmogorov energy spectrum for the velocity. 
The magnetic field is likely to follow in the early kinematic phase Ej^{k) ~ fc 3 / 2 , the Kazantsev solution to the prob- 
lem (see [21] for review). In the case when Sv << Sb, the opposite is true: the structure function cubic in the magnetic 
fluctuations scales like r, in which case again a Kolmogorov spectrum might be plausible, this time for the magnetic 
energy Em, and thus likely for the total energy dominated by Em- Finally, when Sv ~ 5b, we cannot conclude directly 
from these laws the scaling of energy spectra, but we observe that velocity-magnetic field correlations will play an 
essential role, as stressed in |116j (see {VI). In other words one should not factorize the correlators in equations (14 
15), although that would lead of course to K41. 



V. EXPLICIT EXAMPLES OF NON-UNIVERSAL BEHAVIOR IN MHD 



Let us end this section by noting that there are a few known examples of MHD turbulence with different energy 
spectral indices, either in the context of reduced MHD corresponding to a simplification of the MHD equations when 
a strong field is imposed [36], or using the full MHD equations with a strong Bg, in the latter case finding either a 
K41 or an IK spectrum [101] . 

There is also a recent example where all three spectra (K41, IK, WT) are observed when changing in some controlled 
and seemingly simple fashion the initial conditions [77 . Three different runs have been performed, each with the same 
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initial velocity field, specifically the Taylor-Green vortex, a configuration corresponding to laboratory experiments 
between two counter-rotating cylinders, see e.g |97) . Three different initial conditions were taken for the magnetic 
field, with the same magnetic energy and magnetic helicity (zero in that case) and comparable cross-correlations 
(between and 4% in relative terms, i.e. when normalized by the energy). Also equal initial kinetic and magnetic 
energy were taken, and there was no uniform magnetic field. Thus, at t — 0, ta — tjvl, where here ta here is based 
on the rms magnetic field. These computations are performed on equivalent grids of 2048 3 points, imposing the 
four-fold symmetries of the Taylor-Green flow to both the velocity and the magnetic field, taken initially to have 
such symmetries as well. Finally, there is no forcing and the magnetic Prandtl number is equal to unity. What 
distinguishes the runs, which end up with three different spectra, is that, in particular in the case where WT is 
observed, the magnetic energy grows significantly at the expense of the kinetic energy at the largest modes, leading 
to ta(() << tjvl(^) for I ~ L . One should note also that the critical balance hypothesis, even when generalized 
to a non- unity value for the ratio of time scales, is not fulfilled in these simulations [77]. The scale-dependence of 
the time ratio is found to be compatible with the energy spectra scaling obtained for these three runs. In fact, both 
in the computations in [77] and in the parametric study in |141j . the attractive solution seems to be one in which 
quasi-equipartition between kinetic and magnetic modal energies obtain at all wave numbers except in the largest 
scales, and it is these large scales that determine in fine the power-law solution of the dynamics (see Fig. |4j and 

Ea am). 

However, the statistical equilibria for these three initial conditions are the same since the invariants are (almost) 
identical. So, either universality in MHD breaks down, presumably in different sub-classes (possibly K41, IK and 
WT), or else there are hidden invariants of the problem that we are not taking into account (also see [132] for a 
discussion of universality). One ought to consider higher-order correlators such as tetrads involving the Lagrangian 
dynamics of four spatial points [32] , or higher-order correlators involving skewness and kurtosis for example. Is this 
the solution to this conundrum? Or is it that these numerical simulations were done at too low a Reynolds number, 
although the grids were the largest MHD runs performed at the time, with k max /k m i n — 700, where k max , k min 
are the maximum and minimum wavenumber of the computation (corresponding respectively to the size of the box 
and the mesh size)? Since no Bq is imposed, a moderate resolution does lead to a strong large-scale field, but which 
cannot be considered as quasi-uniform still, from the stand-point of the smallest resolved scales. Furthermore, it is 
perhaps important to understand that when one imposes B (or rotation or gravity for that matter), the flow 
is anisotropic at all scales, even if isotropy is recovered at small scales in some cases, like for rotating turbulence; 
whereas, when Bo = 0, there is no preferred direction that can be identified, although one can define a local mean 
Bo,ioc by averaging the magnetic field b in a ball of diameter the integral scale of the flow. Isotropy does recover on 
average at large scale in that second case. 

Another possible reason for these computations to differ in their energy scaling is that they are decaying flows and 
temporal averages cannot be taken, except near the peak of dissipation. Will the forced cases follow the same scaling 
as the decay case, or not? Will quasi-equipartition at the largest scale be achieved on average, or will the dominance 
of magnetic energy in one particular case leading to the weak turbulence regime be able to maintain its status for all 
times? Preliminary results indicate that the non-universality is preserved in the forced case |19j . but these are open 
questions in need of further investigation, the answer to which, though, may depend on the correlation time of the 
forcing as well. 

Finally, note that these exact laws can be used to inverse engineering data and deduce what is the heating rate in 
the Solar Wind due to a turbulent cascade, as done for example in [50], a cascade for which there has been ample 
evidence since the pioneering work done using Voyager data [83] (see also |139| ). 

VI. THE ROLE OF ALIGNMENT BETWEEN THE VELOCITY AND THE MAGNETIC FIELD 

The considerations of the preceding section were done under the (implicit) hypothesis that the correlations between 
the velocity and the magnetic field were weak. However, it has been known for a long time that such correlations 
grow in time (see e.g. |84l I113j ). It was shown in [56) using a simple dimensional argument, that attributing spectral 
indices to the E^ spectra, one could deduce that m + + m~ = 3, with m + = mT =3/2 in the uncorrelated case. 
The reason why both m + and m~ appear in this dimensional analysis is because one Elsasser variable is advected by 
the other one, hence both fields will appear since there is no self-advection. 

This analysis was confirmed in numerical simulations of MHD turbulence both in 2D and in 3D, although in the 
latter case with moderate resolution. Moreover, not only does the global correlation of the flow, which can be viewed 
as the ratio of two invariants, viz. Hq/Et, increases, so does the local correlation which can be defined either as 
p/(x) = v ■ h/[v 2 + b 2 } or pa{x) = cos(v,b), the latter indicating that this is a geometrical property of the flow, 
namely the tendency for the velocity and magnetic field to align or not. This local (as opposed to global) alignment 
was observed in two dimensions [86] as well as more recently in 3D [85], with large plateaux of highly aligned (or 
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anti-aligned) fields and strong gradients in the reconnection regions. It was also shown that other alignments occur, 
between the velocity and the vorticity, as well as between the magnetic field and its current [133] . 

The weakening of cross-correlation in the region of strong gradients was already hypothesized phenomcnologically 
and found numerically using second-order closures [5S] EZ] an d confirmed with 2D-DNS soon afterwards [50] . As 
noted in [155 , it may have important consequences for reconnection since in that domain the interactions can be 
strong (see also |154] for an analysis of Solar Wind data in terms of the anisotropy of the ± Elsiisser variables) . 

The weak turbulence analysis gives an unambiguous answer, namely that m + + mT = 4, which corresponds in 
the closure case to m + + mT = 3 for the isotropic spectra, as reported in the simulations. It also gives the ratio 
of the generalized Kolmogorov constants appearing in front of the E ± spectra as a function of the imposed degree 
of alignment in the flow, with C+ = C_ ~ 0.585 at zero correlation, whereas these two constants ~ 0.19 at high 
correlation [44] . Moreover, the energy fluxes in that latter case differ by more than one order of magnitude. 

However, besides these regimes, the question of what are the asymptotic high Reynolds number spectra in MHD 
turbulence in three dimensions in the presence of strong correlations between the velocity and the magnetic field, 
is highly debated at the present time. Note that MHD in 2D can be viewed as being weak in some sense since it 
is the simplest approximation to the case of 3D MHD in the presence of a strong imposed field. The dynamics to 
lowest order becomes 2D. Or else one can consider the next approximation of the so-called reduced MHD. However, 
in the plane perpendicular to B , interactions can be strong again; and in the absence of a strong imposed field, the 
large-scale isotropy leads to interactions of packets of turbulence in all directions, feeding isotropy to smaller scales 
until possibly the magnetic Taylor scale [90] . 

On the one hand, there is a simple if ad-hoc way to estimate the third-order correlators appearing in Eqs. (Jl4 151, 



as done in [15j : introduce a dimcnsionlcss field, such as an angle, and a natural choice is of course the angle 9 between 
v and b. Thus, the following relation is still dimensionally correct 

(66v 2 5b) ~ r , 

and can be viewed as incorporating the effect of the alignment between velocity and magnetic field through 9. Using 
now the hypothesis of the simplest solution where all variables scale in the same way, one obtains 5v ~ 5b ~ 9 ~ r 1 / 4 , 
compatible with an anisotropic IK spectrum, E(k) ~ kj 3 ^ 2 spectrum, as advocated for some time by Boldyrev 
[15] (see also [25] [52] for another model with similar scaling taking into account the Lagrangian point of view). 
This spectrum is corroborated by moderate resolution DNS in the presence of a strong uniform magnetic field (see 
[8"T] and references therein). Note that, in the forced case, long-time integration can be performed leading to good 
statistics when averaging over the duration of the statistically steady state. For example, it is shown in |107j . using 
the incompressible reduced MHD equations on a grid of 1024 2 x 256 points, a narrow forcing in the perpendicular 
direction and a wide forcing in the parallel direction, with a resulting global correlation coefficient of ~ 0.8, that the 

3/2 

amount of alignment is irrelevant to the final energy spectra which are all found to be proportional to kj_ , with 
"pinning" (joining and possibly crossing) at the dissipation scale, and thus compatible with ideal dynamics (see fjll|. 
Note that another solution, beyond m + + mT = 4, m + ^ to - , is that in fact, at high Reynolds number, the E ± 
symmetry is recovered and m + — mT = 2 for all degrees of alignment between the velocity and the magnetic field, as 
advocated in [8]. 

One of the point that may have to be taken into consideration is based on the importance of the actual values of 
R ± . Indeed, when taking a very strong global degree of alignment, say positive, one has Z$ » Zq ; hence R + could 
be small unless the viscosity itself is extremely small. It is not clear whether these ± Reynolds numbers actually play 
a role in the overall physics, based on the dynamical v and b fields. A similar argument could be made in fluid helical 
turbulence when using the ± variables linked to the eigenmodes of the curl operator. As pointed out in 30 , one may 
be led to the wrong conclusion concerning the dynamics of Hy when considering the Reynolds numbers based on the 
two ± helically polarized waves of the problem. 

Similarly, one question concerns the rates at which the ± energies cascade to small scales, let us call them T±: do 
we have T + ^ T_? Finally, the presence of a bottleneck at small scale at the onset of the dissipative range could be 
an issue as well; this accumulation of energy can be rather prominent for fluids, less so for MHD, a point attributed 
to the greater degree of non-locality of nonlinear interactions in MHD turbulence [27J [3S] ESI HE] ■ 

At this point, it is difficult to conclude and we can consider this problem of strongly aligned MHD turbulence, 
as observed in several instances in the Solar Wind (see e.g. (140] and references therein), as being open from a 
theoretical and numerical point of view. One way forward may be to use models of turbulence to mimic the effect 



of higher Reynolds numbers (see [VIII). Models of small-scale MHD turbulence with non-zero cross helicity have 
been developed |155] and could lead to improvements in the dynamo problem of generation of magnetic fields and in 
reconnection regions as well. 
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VII. CAN THE DYNAMICS OF MAGNETIC HELICITY PLAY A ROLE FOR ENERGY SCALING? 

Kinetic helicity does not seem to play a significant role in fluid turbulence, although the small-scale vortex filaments 
have been known for a long time to be fully helical (see |95j for a review). It may be important in the case of rotating 
flows (see [94J and references therein) , or in the dynamics of the atmosphere, where it is measured in so-called supercell 
storms and hurricanes (see e.g., [US]). The exact law resulting from the flux conservation of kinetic helicity can be 
found in [55j . and that for magnetic helicity was derived in [119 . In terms of the components of the magnetic potential 
and of the electromotive force £ = v X b, it reads: 

{[£ (x) x A(x')]i) = +\e HM r , 



with x' = x + r and £h m — DHm/DT; it can also be written as: 



K(x)£^(x)^(x')) - (6 i (x)E 1 «,(x)^(x')) 



1. 



This gives new constraints on the dynamics of MHD, the consequences of which have not been explored yet. What can 
be the role, then, in the dynamics of MHD turbulence, of magnetic helicity? We saw in fjlljthat its role is essential: it 
is the culprit at the origin of the lack of complete equipartition between kinetic and magnetic energy, and its spectrum 
in the ideal case (when not trivially zero) could be seen as governing all other spatial dynamics except for the kinetic 
energy spectrum. 

Dimensional analysis a la Kolmogorov leads for the magnetic helicity spectrum to a power-law solution HM{k) ~ k~ 2 
|121) . a spectrum observed in the inverse cascade of magnetic helicity in turbulence closures as well as in moderate 
resolution DNS. But in fact other spectra have been observed more recently [791 |9"TI [102] . with Hmik) ~ k~ K 1 with 
k s» 3.3. The spectrum in |121) was obtained under the assumptions of separation of ranges (direct cascade of total 
energy, inverse cascade of Hm), and of locality in Fourier space of nonlinear interactions. But it has been found 
recently that in fact all scales interact for a variety of turbulent fluids [371 GFJ E3j, even though, in a logarithmic 
discretization of scales (as opposed to linear), such may not be the case [3]. And if indeed the magnetic helicity 
spectrum is quite steep, non-local interactions are bound to play a role. Also note that long-time memory effects may 
lead to recurrence of events, as for example in the case of geo-magnetic reversals [TXil ITS] . 

These non-local interactions may render dimensional analysis invalid, and what we need is another way to proceed, 
and perhaps a paradigm shift. The two safe arguments that can be made rely on ideal dynamics described in JET] and 



on the exact laws analyzed in S IV and their extensions to more complex cases. Beyond these two types of results 
and beyond dimensional analysis a la Kolmogorov, there are a variety of arguments that can be put forward, among 
which one that may turn to be valuable (see [79l [9TJ 1110] ) . 

The idea is to postulate that the nonlinear dynamics in MHD, beyond preserving the exact laws, is trying throughout 
the inertial range to achieve a quasi-equipartition between kinetic and magnetic energy and helicity (but not quite 
because of the ideal constraints discussed in £|TTJ) ; in other words, one can expect that there be a proportionality 
between E v (k)/E M (k) and H v (k)/H M (k). This argument can be shown as well to be compatible with the dynamo 
regime [110] . Such a quasi- Alfvenization is observed in the simulations analyzed in |79] , using hyperviscosity at small 
scale. Similarly, in the decaying runs described in [77] using a code that imposed the four-fold symmetries of the 
Taylor-Green flow, quasi- Alfvenization (or QA) is obtained at all scales but the largest (see Fig. 2 in [77]), whereas the 
so-called critical balance postulated in (SHj is not observed (see Fig. 3, op. cit.) . However, one should be careful here: 
this does not mean that critical balance is ruled out; in fact there is evidence for it in some numerical simulations 81J. 
It simply means that it may not be the only solution, and that quasi- Alfvenization may well be a relevant concept as 
well in most of the small-scale inertial range. 



VIII. CAN MODELS HELP? 



Moore's law gives an increase of resolution by a factor two roughly every six years in three dimensions. Hence, 
high Reynolds number turbulence will remain un-attainable in the near future. There are several ways around this 
difficulty, of course. One is to use and impose symmetries of some flows and reach equivalent resolutions that can 
be as much as 8 times a full DNS (see [III [69] for fluids, and [77] for MHD). One can also filter the numerical data 
at small scale, either simply with a higher power of the Laplacian, or else using some approximation to the sub-grid 
stress tensor (see [571 1100] for reviews of Large Eddy Simulations, or LES). 

Another possibility is to use two-point closure models which have been shown to be quite helpful as a source of 
modeling for high Reynolds number turbulence, as for example the Eddy Damped Quasi Normal Markovian (EDQNM) 
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closure. They yield expressions of transport coefficients that depend on time and on the kinetic and magnetic energy 
and helicity spectra. These coefficients can be used as a model of the unresolved sub-grid scale interactions; as such, 
they were recently applied in a variety of conditions, e.g. inhomogeneous flows [UJ, or MHD flows [7], including in the 
study of the dynamo effect 120 , or when helicity is present [155] . In the EDQNM closure in its simplest formulation, 
only two time-scales are considered to damp the high-order correlations, namely the eddy turn-over time and the 
Alfven time. However, when examining carefully the algorithm that leads to a closure formulation, it was shown in 
[7J that another time needs to be considered as well, namely one built on the defect of equipartition between velocity 
and magnetic field, in which case a better model results. Note also that the EDQNM model in the presence of v-b 
alignment has been written [SS] but probably not exploited in the context of modeling of MHD turbulence. This may 
be an avenue to explore. 

Another avenue is to take so-called scalar models, first put forward in MHD in [FT] and expanded further since 
then in several directions (see e.g. |lllj for a recent application to the dynamo in a rotating fluid at small magnetic 
Prandtl number). 

Finally, let us mention another filtering methodology which may be quite promising in MHD |109| : in that approach, 
the invariants are maintained but in a different norm, Hi instead of L 2 (see e.g. |64j ) . In other words, within the 
model, it is Ejji ~ (l u | 2 + o 2 |oj| 2 ) which is conserved, with a an open parameter of the model taken as the scale 
on which Lagrangian trajectories are averaged. This allows constraining the growth of vorticity and current density, 
and thus reaching higher Reynolds number at a given resolution than the DNS would allow, by up to a factor 6 or 
perhaps more [110] , For example, one finds in |110j a clear indication, for an equivalent resolution of ~ 6000 3 points, 
that quasi Alfvenization is the preferred mode for this given set of initial conditions: the energy ratio EM{k)/Ey(k) 
remains remarkably constant (and rs 2) throughout the inertial range, and with a ratio of Alfvenic to turn-over time- 
scales compatible with the scaling of the energy spectrum. This clearly is worth exploring further, at higher Reynolds 
number as well as in the forced case. 

A combination of all these modeling techniques may prove a valuable tool for studying MHD turbulence at Reynolds 
numbers as encountered in geophysics and astrophysics, in order to explore some of the issues mentioned in the 
preceding Sections. 



IX. CONCLUSIONS 



A lot more needs to be done, using every tool we can, to help unravel what is happening in MHD turbulence. Data, 
both numerical and observational, seems to indicate that, at orders that are higher than the energy spectra, there is 
a notable difference between fluid and MHD turbulence, due to the dynamics of small scales, vortex filaments in one 
case and vorticity and current sheets in the other, which do eventually roll-up as observed in the magnetosphere (see 
e.g. [HI])- Whereas it is difficult to measure small discrepancies in scaling exponents, it might be possible to explore 
the physical differences that lead to critical balance versus quasi- Alfvenization, as discussed in §VII| in the context of 
magnetic helicity which may well play a determining role, possibly as important as the degree of alignment between 
the velocity and the magnetic field. Perhaps, a study of Lagrangian dynamics as performed for fluid turbulence and 
recently in MHD as well (see "2"5J US]) will help understand the detailed properties of small-scale in MHD. 



In fact, when discussing the exact laws written in Eqs. (14 15), more possibilities could take place beyond the 



relative ordering of velocity and magnetic fluctuations, in introducing the degree of alignment between v and b on 
the one hand, and between A and b on the other hand. According to which alignment is greater, one regime can 
dominate over another one. Also worthy of noting is the fact that, in the ideal dynamics of MHD turbulence, the 
cross-correlation and magnetic helicity spectra are proportional with He ~ Ic 2 Hm, with a coefficient of proportionality 
depending only on the three Lagrange multipliers associated with the three invariants. Since one has postulated and 
observed an inverse cascade of magnetic helicity to the large scales (see e.g. [4*21 11211 1123| ). then would it be possible 
to envisage that, for some ratios of the a,f3,j multipliers, the cross-correlation could also grow at large scale (even 
though, in that case, fc — > 0), implying a high degree of alignment of the velocity and magnetic field in the large scales 
of the system? There is some evidence for this in the Solar Wind |140| . but more data needs to be analyzed in detail 
before we can conclude on this point which may be mature for numerical explorations as well (see also equation) 10 )). 
Also, in the case of the interstellar medium, perhaps, with ALMA for example, more detailed information about the 
scaling behavior of interstellar turbulence will be available (see e.g. [40]). 

On the other hand, structures (such as horse-shoe vortices in channel flows, or convective plumes) arise from the 
boundaries, and one may wonder what is the role of such boundaries in determining scaling and statistical properties 
of turbulent flows in general. This question can be extended of course to homogeneous turbulence which naturally 
develops internal boundary and shear layers, and to the effect of the roughness of such boundaries [54j : do current 
sheets (and shear layers) provide a similar corrugation of MHD turbulence and lead to specific scaling laws? 

Numerically, one of the problem is the lack of adequate resolution in three space dimensions: how high do we need 
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to push the Reynolds number to be convinced that we have an asymptotic solution? Statistical mechanics does tell 
us that diverse spectra can emerge. Will K41 still emerge at very high Reynolds? Or does the intuition based on the 
examination of the exact laws stemming from the conservation properties of MHD indicate indeed that three different 
regimes can occur, as seen in some computations? 
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